function caxout = PlotSurf(ODforplot,caxin)

surf(ODforplot,'Marker','.') 
if nargin == 2
    caxis(caxin)
end
caxout = caxis;
shading flat
grid off
axis off
view(0,90);
end